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Manipulation and Generation of Supercurrent in Out-of-Equilibrium Josephson 

Tunnel Nanojunctions 
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We demonstrate experimentally manipulation of supercurrent in Al-AlOaj-Ti Josephson tunnel 
junctions by injecting quasiparticles in a Ti island from two additional tunnel-coupled Al super- 
conducting reservoirs. Both supercurrent enhancement and quenching with respect to equilibrium 
are achieved. We demonstrate cooling of the Ti line by quasiparticle injection from the normal 
state deep into the superconducting phase. A model based on heat transport and non-monotonic 
current-voltage characteristic of a Josephson junction satisfactorily accounts for our findings. 

PACS numbers: 74.50.+r,85.25.Cp,73.23.-b,74.78.Na 
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Nonequilibrium dynamics in superconducting nanocir- 
cuits is currently in the focus of an intense experimental 
and theoretical effort [l|, |2j. In this context, the con- 
trol of the Josephson current in superconductor-normal 
metal-superconductor (SNS) weak links is receiving much 
attention. In these systems supercurrent is manipulated 
by modifying the quasiparticle energy distribution in the 
N region via current injection from external terminals 
Pi- y, [5(- There have been some successful demonstra- 
tions of such out-of-equilibrium SNS junctions [a, Lfl, 13, l2( ■ 
On the other hand, it was predicted lfj, [llj that super- 
current can be controlled in all-superconducting tunnel 
structures as well. In this case quasiparticle injection 
12, |l3|, [lj, ll5j, LLa [l7| , leads to intriguing features pecu- 
liar to out-of-equilibrium superconductors. 

In this Letter we report on control of the Josephson 
coupling in a small S island by injecting quasiparticles 
from tunnel-coupled superconducting leads. Both su- 
percurrent enhancement and suppression with respect to 
equilibrium, as well as generation at temperatures above 
the island critical temperature were achieved by chang- 
ing the quasiparticle injection rate. Our findings are 
explained within a model relating the superconducting 
state of the island to the heat flux driven through it upon 
injection. 

Figure [1] shows a scanning electron micrograph of a 
typical structure along with a scheme of the measurement 
setup. The core of the sample consists of a SIS IS con- 
trol line, i.e., a titanium (Ti) superconducting island (S ) 
symmetrically connected at its ends via AlOa; barriers (I) 
with normal-state resistance i?x each to two aluminum 
(Al) superconducting reservoirs (S). Two additional Al- 
AlO^-Ti probe junctions, with normal-state resistance 
i?j each and placed in the center of the island, are used 
to measure the Josephson current (Ij) of their nomi- 
nally symmetric series connection. Our samples were 
fabricated by electron beam lithography and two-angle 
shadow-mask evaporation. The measurements were per- 
formed in a dilution refrigerator at sub-kelvin tempera- 
tures measured with a Ru02 resistor calibrated against 




FIG. 1: (color online) A typical structure (Sample B) show- 
ing a schematic of the measurement setup. In the middle, a 
Ti superconducting island (S ) is connected to four Al elec- 
trodes (S) through tunnel junctions. Ij denotes the Joseph- 
son current flowing through the two inner Al-AlO^-Ti tunnel 
junctions. 



Coulomb blockade thermometer [2| . The experiment con- 
sists of measuring at different bath temperatures (Xbath) 
the current-voltage characteristic (I pr vs V pI ) of the se- 
ries connection of the central SIS Josephson junctions 
while imposing a fixed voltage (V; n j) across the lateral Al 
reservoirs. As we shall show, this will lead to a change in 
temperature of S which determines the dynamics of the 
Josephson junctions. 

Figure [2] shows the electrical characterization of two 
structures, in the following referred to as Sample A 
(whose essential parameters are i?x — 1.43 kfi, Rj ~ 2.8 
kfi and a 45-nm-thick Ti island of area 250 x 2550 nm 2 ), 
and Sample B (with R T ~ 710 0, R 3 ~ 1.56 kfi and 
a 40-nm-thick Ti island of area 650 x 1500 nm 2 ). The 
critical temperature (T c ) of the Ti layer is ~ 500 mK for 
Sample A and ^210 mK for Sample B. Panels (a) and 
(b) display the low-temperature I pl - vs V pr characteris- 
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FIG. 2: (color online) I pr vs V pl for various values of Vinj: (a) 
Sample A, T bath = 62 mK; (b) Sample B, T bath = 35 mK; (c) 
Sample A, T bath = 510 mK; (d) Sample B, T bath = 250 mK. 
The curves are vertically offset for clarity. 



tics of Sample A and B, respectively, for several values 
of the injection voltage Vinj- Each characteristic corre- 
sponds to a different Vinj, and the curves are vertically 
offset for clarity. In equilibrium, at Vinj = 0, the super- 
current manifests itself as a peak around zero bias in the 
current- voltage characteristic. As will be explained with 
further details, upon increasing the injection voltage the 
supercurrent behaves non-monotonically, being initially 
suppressed, then showing typically two peaks at inter- 
mediate injection voltages. Further increase of Vnj leads 
to a monotonic supercurrent decay, and to a complete 
quenching for Vnj > 400 /xV. In Sample B the peak ampli- 
tude is enhanced by almost a factor of three with respect 
to equilibrium. The supercurrent response in the high- 
temperature regime [see panels (c) and (d) for Sample A 
and B, respectively] is different. In particular, while the 
equilibrium supercurrent is already vanishing, as Tb a th 
exceeds the critical temperature of the Ti island, it is 
generated by increasing Vi n j at an injection voltage which 
cools S from the normal into the superconducting state. 
This occurs thanks to hot quasiparticlc extraction pro- 
vided by the Al reservoirs [2j, [lfj, [ll| . By increasing V; n j 
even further leads to another maximum of supercurrent 
followed by full suppression. Although somewhat differ- 
ent in terms of characteristic parameters, both samples 
show similar behavior. 

The full dependence of the maximum supercurrent 
I max on Vi n j at different bath temperatures is displayed 
on the left axis of Fig. [3ja) and (b) for Sample A and 
B, respectively. I max is defined as the average between 
the amplitudes of positive and negative peaks of I pi . It 
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FIG. 3: (color online) (a) Left axis: I max vs Vinj at three dif- 
ferent T ba th for Sample A. Right axis: Injector characteristics 
Jinj vs Vi n j at the same bath temperatures, (b) The same as 
in (a) for Sample B. (c) 7 max vs T ba th for three different values 
of Vinj for Sample A. (d) The same as in (c) for Sample B. 



is a symmetric function of Vinj based on electron-hole 
symmetry, so that just the dependence on positive Vj n j 
is shown. As we shall show in the following, the features 
present at small Vbi as in the supercurrent response are 
related to the shape of the current-voltage characteris- 
tic of the injectors (Jj n j vs Vinj), shown on the right axis 
of Fig. [3]^a) and (b) for the same Tbath- In particular, 
in addition to a current enhancement around Vinj = 
originating from Josephson coupling in the lateral SIS 
junctions, the curves at lower Tbath show a marked peak 
centered around the middle of the characteristic which 
disappears as soon as S undergoes a transition into the 
normal state. 

Figures [3Jc) and (d) show the J ma x vs Tbath character- 
istic for Sample A and B, respectively, at three different 
values of L m j . For Vj n j = (open squares) the equilibrium 
supercurrent saturates at low Tbath where it obtains val- 
ues as high as ~ 2.35 nA and ~ 4.8 nA for Sample A 
and B, respectively, while it is gradually reduced by in- 
creasing the temperature, being completely suppressed 
at Tbath — 500 mK and ~ 210 mK, i.e., at the critical 
temperature of Sample A and B, respectively. The low- 
temperature supercurrent amplitudes are suppressed in 
both samples by about an order of magnitude as com- 
pared to the Ambegaokar-Baratoff theoretical prediction 
[181 ] . This is however expected for ultrasmall Joseph- 
son tunnel junctions influenced by environment fluctu- 
ations [19|, [2fj|. For a chosen injection voltage, e.g., at 
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325 fj,V, / max saturates at low Tb a th at ~ 2.55 



nA for Sample A, and ~ 10.8 nA for Sample B. The 
maximum supercurrent survives under injection up to 
Tbath — 630 mK for Sample A and ~ 360 mK for Sample 
B, i.e., well above the equilibrium critical temperature. 
This means that we can cool the samples by quasiparticle 
current from the normal into the superconducting state. 
Also shown is the temperature dependence at the opti- 
mized bias voltage (Vopt) which maximizes / max (solid 
dots). 

Our observations of non-monotonic dependence of the 
probe supercurrent on bias voltage and of a peak in 
the current in the middle of the superconducting gap 
can both be explained qualitatively within a very sim- 
ple model. The key observation is that the bias voltage 
dependence of the current of a single injector junction 
is non-monotonic because it can be carried by Cooper 
pairs (supercurrent around zero voltage) and by quasi- 
particles (near and above the gap voltage). Then, as a 
function of bias voltage Vinj across the two injecting junc- 
tions, the evolution of the voltage across each individual 
junction is as follows [see the energy-band diagrams in 
Fig. UK a)]. At around zero bias, both junctions carry 
supercurrent, seen as a peak in the current-voltage char- 
acteristic. Thereabove one of the junctions, i.e., the one 
with smaller critical current [for instance, left (L) injec- 
tor in Fig. Ufa)], switches into the quasiparticle branch: 
the total voltage then equals that across this "weaker" 
junction, while the other one remains in the approxi- 
mately zero voltage supercurrent branch. In this situ- 
ation, when the voltage is approximately (As — Ag/)/e 
|2j, llCj, llj| , there is increase of Josephson critical current 
of the probe junctions, thanks to enhanced cooling power 
(<2l) due to quasiparticle current in the L junction, i.e., 
Ql 7^ 0. Here A s g > is the BCS energy gap in S (S ). 
In the middle of the gap region at (As + As')/e, one 
of the junctions reaches the steep onset of quasiparticle 
current leading to a peak in control current. Above this 
bias, also the second junction [i.e., right (R) injector in 
Fig. HJa)] switches into quasiparticle branch providing 
finite cooling power, i.e., Qr ^ 0. Now the voltage is 
divided approximately equally across the two junctions, 
and at intermediate voltages above (As + As')/e cool- 
in g po wer is small until it maximizes at 2(As — As')/e 
|2l.ll0llll| resulting in another maximum in probe super- 



current. The final increase of current ii n j in the control 
junctions occurs at 2(As + Ag')/e, where both junctions 
have an approximately equal voltage corresponding to 
the onset of quasiparticle current: this results in large 
current, heating of the S island, and subsequent quench 
of the probe supercurrent. 

A more quantitative analysis can be carried out as fol- 
lows. The total electric current flowing through left and 
right injectors can be written as I.\ = Ij ' + /qp R , 

is the Ambegaokar- 



(a) S(a s ) 

V 12 



Q L * O r = 



' L S'(A S .) R 



S(A S ) 



-VJ2 



(fc£V| nJ SlA s +A s .)/e 



S(A S ) 



VJ2 



V m > (As+AsO/e 



■(b) 



/'' Tj£"250mK f "_ 



•Lr 



200 mK 



130 mK 



35 mK 




' 71 .(mK) 200 

i bath* ' 

1 30 

I , 35 



200 
K>V) 



400 



200 400 



FIG. 4: (color online) (a) Energy-band diagram of the biased 
SIS IS control line. Also shown on the right side is the volt- 
age drop Vh across left interface, (b) Left axis: 7j vs Vi n j 
calculated for a few values of Tb a th- Right axis: Ji n j vs VS n j 
calculated for the same Tbath- The curves are vertically offset 
for clarity, (c) Calculated T a vs Vi n j for the same Tb at h as in 
(b). The horizontal line indicates Ti critical temperature of 
Sample B. 



Baratoff critical current of the injectors [l8j, Vrr is 
the voltage drop across L(R) interface [see Fig. 0Ja)], 
while 7 q L ' R = ±^ 7 /AAAs(eL,R)AA S '(e)[/o(?L,R,T bath ) - 

/o(e, T c )] is the quasiparticle current. Here, cl,r = 
e T e^nj/2, e = e - e(V inj /2 - Vl), / (e, T) is the Fermi- 
Dirac function at temperature T, and Af s s > (e) is the 

smeared density of states of S(S ). In particular we 
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where Jj ' ^ for Vl,r 



set Af s g / (e) = |Re[(e 

where r g g / accounts for quasiparticle states within the 

gapinS(S') [2]. 

The voltage drop across L(R) interface resulting from 
biasing with V4 n j follows from the conservation of the to- 
tal current, i.e., 7 mj = 7 mj with V L + Vr = V in j. The 
solution for Vl [shown on the right side of Fig. Ufa)] is 
2V L /V hli - 1 = ±1 for < V inj < (A s + A g »)/e, and 
Vl = Mnj/2 for V m j > (As + A g ')/e, meaning that only 
one junction is initially in the dissipative regime [L (R) 
junction in the upper (lower) branch]. The threshold for 
equal voltage division is Vinj — (As + A s ')/e and it de- 
pends only marginally on the asymmetry between the 
two injector junctions. 



The supercurrent of the probe junctions depends on 
the quasiparticle distribution in S under voltage bi- 
asing [lfj|, LUJ- Strong electron-electron interaction 
drives the electron system in S into local thermal 
(quasi)equilibrium described by a Fermi-Dirac function 
at an electron temperature T which may differ from 
Tbath 2]. The maximum Josephson current flowing 
through the central SIS junctions is given by l(llll| 
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above expressions we set As = As (Tbath) and A g ' — 
A g '(T c ). Equation (f) shows that Ij is controlled by T 
once Tbath is fixed. Under bias voltage Vinj the heat cur- 
rent (Ql.r) flowing from S to S through L or R interface 
is given by [Tfl, 12JJ 



f 



?l,r = ^- / deeN s (£ L ,R)Af s > {e) 
x[/o(e, T c ) - / (e L ,R, Tbath)]- 



(2) 



at Vi n j — 150 /xV. Further increase of bias voltage leads 
again initially to heating, then cooling, and eventually 
heating above T c for large Vi n j . At the bath temperature 
of 250 mK, T c starts to decrease monotonically, initially 
driving S into the superconducting state, and showing 
the same behavior as at lower Tbath- 

In conclusion, control of Josephson current as well as 
its generation at bath temperatures above the critical 
one were achieved by varying quasiparticle injection into 
a small superconducting island. Our results are success- 
fully described within a model relating the superconduct- 
ing state of the island to the heat flux originating from 
quasiparticle injection. From the practical point of view, 
our experiment demonstrates that quasiparticle injection 
can cool a metal wire from its normal state deep into the 
superconducting phase. 

We acknowledge financial support from the EU 
Large Scale Installation Program ULTI-3 and from the 
NanoSciERA "NanoFridge" project. 



T c is then determined by solving the energy-balance 
equation Q h {V inj , T ba th, T^)+Q R (V inh T batn , T e ) = 0. We 
neglect the electron-phonon interaction contribution in 
the energy-balance equation which would lead to small 
corrections only. The probe supercurrent is then deter- 
mined by the electron temperature T e established in S 
by biasing the control line. 

For comparison with the experiment we chose the given 
parameters of Sample B, T c = 1.2 K and depairing pa- 
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The injector current- 
voltage characteristics calculated at different Tbath are 
displayed on the right axis of Fig. [4|b). In addition to 
Josephson coupling vanishing at Tbath > T c , the current 
shows a peak centered in the middle of the characteris- 
tic, as observed in the experiment [see right axis of Figs. 
EJa-b)]. The 7j vs Vj n j characteristics are displayed on 
the left axis of Fig. [4jb) for the same Tbath values. The 
supercurrent curves of Fig. [3ja-b) resemble those of the 
model presented in Fig. HJa), apart from details that we 
attribute to the oversimplified thermal model. 

Figure HJc) shows the electron temperature T e cal- 
culated from the energy-balance equation for the cor- 
responding bath temperatures. For Tbath < 200 mK 
the electron gas is initially heated, inducing supercurrent 
suppression at small bias voltages. Such heating stems 
from subgap current in a tunnel junction [2j, |22j, |23| . By 
increasing V; n j further, the electron temperature starts to 
decrease, thanks to quasiparticle cooling [2], ]2l| provided 
by the larger gap superconductor (S), and is minimized 
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